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Abstract 

The issues of quintessence and cosmic acceleration can be discussed in the 
framework of theories which do not include necessarily scalar fields. It is possi- 
ble to define pressure and energy density for new components considering effective 
theories derived from fundamental physics like the extended theories of gravity 
or simply generalizing the state equation of matter. Exact accelerated expanding 
solutions can be achieved in several schemes: either in models containing higher 
order curvature and torsion terms or in models where the state equation of matter 
is corrected by a second order Van der Waals terms. In this review, we present 
such new approaches and compare them with observations. 
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1 Introduction 



One of the recent astonishing result in cosmology is the fact that the universe is acceler- 
ating instead of decelerating along the scheme of standard Friedmann models as everyone 
has learned in textbooks [1, 2, 3]. Type la supernovae (SNe la) allow to determine cos- 
mological parameters probing the today values of Hubble constant Hq in relation to the 
luminosity distance deduced from these stars used as standard candles [4, 5]. 

For the red-shift z < 1, by the luminosity distance cIl ~ Hq^z + (1 — q )z 2 /2], the 
observations indicate that the deceleration parameter is 

-i < go < o (l) 

which is a clear indication for the acceleration. 

Besides, data coming from clusters of galaxies at low red shift (including the mass-to- 
light methods, baryon fraction and abundance evolution) [6], and data coming from the 
CMBR (Cosmic Microwave Background Radiation) investigation (e.g. BOOMERanG)[7] 
give observational constraints from which we deduce the picture of a spatially flat, low 
density universe dominated by some kind of non-clustered dark energy. Such an energy, 
which is supposed to have dynamics, should be the origin of cosmic acceleration. 

If we refer to the density parameter, all these observations seem to indicate: 

fi Af ~0.3, tt A ~0.7, fi fc ~0.0, (2) 

where Qm is the amount of both non-relativistic baryonic and non-baryonic (dark) mat- 
ter, A is the dark energy (cosmo logical constant, quintessence,..), k is the curvature 
parameter of Friedmann- Robertson- Walker (FRW) metric of the form 



ds 2 = dt 2 -a(tf 



<I,J +r 2 dQ 2 



1 — kr 2 



(3) 



here a(t) is the scale factor of the universe 1 . 

The relations (2) come from the Friedmann-Einstein equation 



87rG k 

* 2 = —p--i (5) 



where H — — is the Hubble parameter. Specifically we have 
a 



_ 8nG N pM n _ 8nG N p A A k 
Um ~ 3H 2 ' ^ A "^i^"3i^' k ~ J^H) 2 {b) 



1 The deceleration parameter can be given in terms of densities parameters. We have for FRW models 

9o = ~ = ^(3 7 +1)Om-^a, (4) 
where 7 is the constant of perfect fluid state equation p = 7/?. 
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for the various components of cosmic fluid. Dividing by H 2 , we get the simple relation 

i = n M + + fife (7) 

and then, through observations Eq. (2). 

After this discovery, a wide debate on the interpretation of data has been developed, 
in particular about the explanation of this new A component, featured as a negative 
pressure fluid, needed for the best fit of observations. It is clear that one of the most 
important challenge for the current and future researches in cosmology, beside the tradi- 
tional issues as large scale structure, initial conditions, shortcomings of standard model 
and so on, is understanding the true fundamental nature of such a dark energy. 
Several approaches has been pursued to this discussion until now. Neglecting some in- 
teresting but more philosophical speculations (i.e. so called Anthropic Principle see [8]), 
the most interesting schemes can be summarized into three great families: cosmological 
constant, time varying cosmological constant and quintessence models. 
The first approach is related to a time-independent, spatially homogeneous component 
which is physically equivalent to the zero-point energy of fields [9, 10]. Such a scheme is 
of fundamental importance since fixing the value of A should provide the vacuum energy 
of gravitational field. 

Deriving such value is important also in the framework of the so called No-Hair conjecture 
whose issue is to predict what will be the fate of the whole universe as soon as, during 
the evolution, A will remain the only contributions to the cosmic energy density 2 [14, 15]. 
This consideration implies shortcomings in such an approach. The predicted value of A 
parameter is very tiny compared to typical values of zero-point energy coming from par- 
ticle physics (at least 120 orders of difference), this discrepancy gives rise to the so called 
problem of cosmological constant. Unfortunately, if we were able to explain the small 
value of the cosmological constant, it would not be enough. A second puzzle, likewise 
not easily explainable, is the why matter density and cosmological constant component 
have today comparable values (coincidence problem). In order to solve these problems, 
many authors have investigated others alternative schemes. 

For example, a time-varying cosmological constant has been introduced, in particular in 
relation to the inflationary paradigm. 

The main feature of inflationary approaches is the development of a de Sitter (quasi-de 
Sitter or power law) expansion of early universe. 

To connect such an expansion to structure formation and to a forthcoming de Sitter 
epoch (as requested by No-Hair conjecture) [18, 19] we need a cosmological constant 
which value acquired a great value in early epoch, underwent a phase transition with a 
graceful exit and will to result in a small remnant to ward the future [20]. This scheme 

2 A more precise formulation of such a conjecture is possible for a restricted class of cosmological 
models, as discussed in [14]. We have to note that the conjecture holds when every ordinary matter 
field, satisfies the dominant and strong energy conditions [16]. However it is possible to find models 
which explicitly violate such conditions but satisfies No-Hair theorem requests. This conjecture can be 
extended to the case of time-varying cosmological constant [17] 
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provides a mechanism to overcome the coincidence problem. 

In fact by a dynamical component, it could be more natural to achieve comparable values 
of energy density between cosmological fluids at a given epoch. 

The straightforward generalization is to consider an inhomogeneous and un-clustered 
cosmological component. From this view point we arrive to the definition of quintessence 
[11], which is a time- varying spatially inhomogeneous fluid with negative pressure. 
It is clear that these three approaches are strictly linked each other. Constant vacuum 
energy density, quintessence and time varying cosmological constant converge to the 
limit 7 — > — 1, since inhomogeneities are not revealed for dark energy component at the 
observable scales. All the authors assume quintessence as just a time- varying cosmolog- 
ical component. In summary all the above schemes claim for an ingredient which is a 
non-clustered energy counterpart which we need to explain observations (in particular 
acceleration of cosmic fluid). 

In order to solve the cosmological constant and coincidence problem, people widely take 
into account scalar fields which are a form of matter capable of giving rise to a state 
equation with negative pressure. The key ingredient of such a dynamics is the form 
of interaction potential which several times is phenomenological and unnatural since it 
cannot be related to some fundamental effective theory. 

In this review, we want to show how quintessential scheme can be achieved also without 
considering scalar fields as usually discussed in literature. 

First, we will show how an accelerated behaviour of cosmic fluid can be achieved assum- 
ing a more physically motivated equation of state for matter: in particular taking into 
account second order terms of Van der Waals form. This approach seems very intriguing 
in relation to the possibility of mimicking some effects, as phase transitions, which char- 
acterize standard matter behaviour. 

From another point of view, we investigate the extended theories of gravity, that is the- 
ories of gravity which generalize the standard Hilbert-Einstein scheme. In this setting, 
it is possible to obtain a cosmological component by a geometrical approach, relating 
its origin either to effective terms of quantum gravity (higher order theories of gravity) 
or to, for example, effective contributions of matter-spin interaction which modify the 
space-time geometry (torsion). 

The paper is organized as follows. Sect. 2 is devoted to a review of principal results 
coming from observational investigations as supernovae, cosmic microwave background 
and Sunyaev-Zeldovich method for clusters of galaxies. The most peculiar evidence com- 
ing from the bulk of data is the fact that today observed universe is accelerating. 
In Sect. 3, we outline the scheme of quintessence in the framework of scalar field approach. 
Sect. 4 is devoted to the Van der Waals quintessence. It is interesting to note that the 
only request of taking into account a more physically motivated equation of state gives 
rise to the acceleration, matching the observations. 

In Sects. 5 and 6, we take into account curvature and torsion quintessence considering 
the presence of geometrical terms into the Hilbert-Einstein action. Such ingredients are 
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essential every time one takes into account effective quantum corrections to the gravita- 
tional field [21]. 

Sec. 7 is devoted to the discussion of results, conclusions and further perspectives of our 
approach. 

2 What the observations really say 

In this section we give a summary of the present status of observations which claim for 
the presence of some form of cosmological dark energy. However the picture which we 
present is far from being exhaustive but it gives an indication of the problem from several 
point of view. 



2.1 The data from SNela 

The most impressive result of current observational cosmology has been obtained by 
the study of supernovae of type la (SNela). The phase of supernova is the product of 
collapse of super-massive stars which degenerate in a very powerful explosion capable of 
increasing the magnitude of the star to values similar to those of the host galaxy. The 
physics of this late stellar phase is still matter of debated [22], however it is possible to 
classify supernovae in some types in relation to the spectral emission. Type la indicates 
supernovae originated by white dwarf carbon and oxygen-rich degenerated through the 
interaction in binary systems. These supernovae are interesting from a cosmological point 
of view since it has been found a characteristic phenomenological relation between the 
amplitude of the light-curve and the maximum of luminosity (Phillips — relation). In this 
way, they can be considered as good standard candles. This feature has an immediate 
cosmological relevance: such kind of supernovae can be detected at high red-shifts (e.g. 
till z ~ 1). 

In order to see how supernovae work as standard candles we have to give some notions 
related to the luminosity distance [2, 23]. 

It is possible to define the luminosity distance di of an astrophysical object as a function 
of cosmological parameters. In tha case of a FRW, spatially flat, metric we have: 

d L (z) = (i + z) f dz'[n M (i + z'f + n K }~ 1 / 2 . (8) 

Jo 

Now, using the well known magnitude-red-shift relation 

fi(z) = 5log-^-d L (z) + 25, (9) 

-no 

we obtain a measure of the distance modulus fi in terms of cosmological parameters. 
Distance modulus can be obtained from the observations of SNela. In fact, the apparent 
magnitude m is measured, while the absolute magnitude M may be deduced from the 
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intrinsic properties of these stars and some adjustment of the Phillips relation as the 
Multi-Color Light-curve Shape (MLCS) method [24]. The distance modulus is simply 
fi = M — m. Finally the red-shift z of the supernova can be determined accurately 
from the host galaxy spectrum. At this point, observing a certain sample of supernovae 
[4, 5] it has been possible to fit different cosmological models which give rise to different 
luminosity distances. 

The two groups SCP (Supernovae Cosmology Project) [4] and HZT (High z Search Team) 
[5] found that distant supernovae are significantly dimmer (nearly half a magnitude), 
with respect to a sample of nearby supernovae, than it would expected in a cosmological 
model with Q M = 1, (that is the standard cold dark matter model). 
Referring to the best fit values, the SCP group suggests a universe with 

0.8Q M -0.6fi A ~ -0.2 ±0.1, (10) 

which gives for a flat model with 

0.28, Vl A ~ 0.72. (11) 

On the other side, HZT constraints are compatible with the ones of SCP; in fact also in 
this case the best-fit value for the flat case is VLm — 0.28. 

From these estimates of matter density parameter, it results that a spatially flat universe 
should be filled by a dark cosmological component with a negative pressure and density 
parameter ~ 0.7. 

As consequence it is easy to obtain a negative value for the deceleration parameter dis- 
cussed above. 

Other results coming from these observations are a best fit value for H ~ 65 km s ec~ 1 Mpc~ 
and an estimate for the age of universe of about 15Gyr. 

2.2 Cosmic microwave background radiation 

Another class of astrophysical observations with great relevance in cosmology is the study 
of CMBR. This radiation has been produced at recombination, when matter is become 
transparent to light and because cosmic expansion has cooled the universe at a today 
temperature of about To = 2.73K. The CMBR can be assumed coming from a far 
spherical shell around us, the so called last scattering surface and can be considered, in 
a first approximation, homogeneous and isotropic. 

Now, from theories of structure formation we know that clustered matter is the result 
of evolution of primordial perturbations acting as seeds. So, it seems obvious to think 
to an influence of such perturbations on the background radiation. We expect primary 
and secondary anisotropies. The first ones should be produced just at recombination as 
an imprint of inhomogeneities on the last scattering surface (Sachs- Wolfe effect, intrinsic 
adiabatic fluctuations, etc..) [1, 2]. The secondary ones are, instead due to the effect of 
scattering processes along the line of sight between the surface of last scattering and the 
observer. 
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By considering only primary fluctuations, we can deal with the matter of early universe 
as a fluid of photons and baryons. We have a competition of gravity and radiation 
pressure effects in this fluid which implies the setting up of acoustic oscillations. At the 
decoupling, these oscillations have been frozen out in the CMBR and, today, we detect 
them as temperature anisotropies in the observed sky. The relevance of such studies is 
related to the possibility of linking the amplitude of fluctuations, and in particular their 
power spectrum, to the spatial geometry of the universe. It can be shown that angular 
scale or, equivalently, the multipole momentum I of the first acoustic peak of spectrum 
can be expressed in terms of spatial curvature. It is possible to show that the relation 
[25] 

200 

l ~v^m (12) 

holds. It gives the first multipole as a function of energy density of spatial curvature. 
In 1992, the COBE satellite showed for the first time that such fluctuations in the CMBR 
are really existing. Starting from such result, others surveys have been performed with 
the aim of probing thermal fluctuations in the sky with a best sensibility. These experi- 
ments, COBE-DMR (COsmic Background Explorer- Differential Microwave Radiometer) 
[26], BOOMERanG {Balloon Observations Of Millimetric Extragalactic Radiation and 
Geophysics)[7], MAXIMA (Millimeter Anisotropy experiment IMaging array )[25], have 
indicated for I the value w 197 ± 6 (BOOMERanG) and w 220 (MAXIMA). The esti- 
mate of spatial curvature density parameter obtained combining the different results is 
ps 0.11 ± 0.07 [27]. It must be stressed that such a number has been obtained taking 
into account a value of Hubble parameter which is deduced thanks to the supenovae 
observations. 

Previous results for Qk were slightly different from zero, but it is possible to reduce them 
near to zero with a reasonable agreement to the others cosmological parameters' esti- 
mates [27]. In relation to this outcome, one can say that the best-fit CMBR results give 
a picture of the universe that is a spatially flat manifold. Such a result is in agreement 
with the case of SNela surveys. 



2.3 Others observational approaches 

In addition to the investigations of CMBR and SNela surveys several astrophysical ob- 
servations can have cosmological relevance, in order to probe structure and dynamics of 
cosmic fluid. 

Among these ones, weak gravitational lensing is a very active research area which is giv- 
ing extremely interesting result [28]. 

The relevance of lensing as a cosmological probe is linked to the fact that different cos- 
mological models give different probabilities for the light coming from very distant object 
to meet gravitational lenses on their walk. This fact can be used to test models with 
observations, in fact, as seen for the supernovae case, the luminosity distance depends 
on the cosmological parameters as the cosmic volume at a specific red-shift. As a result, 
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the counting of apparent density of observed objects, whose actual value is assumed as 
known, provides a test to probe the universe components [29]. Considering a flat uni- 
verse, the probability of a source, at red-shift z, to be lensed, relative to the case of an 
Einstein-de Sitter universe (Qm = 1,^a = 0) is: 

1 Fi c a * 

P lens = ^[l-(l + zT 1/2 ]- 3 l da 

where d A — (1 + z)~ 2 d,L is the angular distance and a* = (1 + corresponds to the 
redshift z*. It must be underlined that such observational approach is frustrated by sev- 
eral uncertainties due to evolution, extinction and so on. Taking into account such bias 
sources, it is possible to obtain an upper limit for cosmological energy density component 
Q A which is < 0.7 [30, 31, 32]. 

We can cite also other recent studies conducted by weak lensing approach: in these cases 
the amplification is dependent on cosmological models and the data indicate always a 
nonzero A — term [33]. 

Other important tests for cosmological parameters are related the evaluation of Qm- 
Studies on Qm have provided values ranging between 0.1 and 0.4. A result drastically 
larger than the density parameter for the baryon matter as inferred from nucleosynthesis 
Q bar w 0.04 [34, 35]. 

Typically the estimates of Qm are performed weighting the mass of galaxy clusters in 
relation to their luminosity and extrapolating the results to the whole universe. Adopt- 
ing this strategy, and using the virial theorem for the dynamics of clusters, a value of 
Qm = 0.2 ±0.1 is obtained [36]. Also others collaborations, which weight clusters using 
gravitational lensing of background galaxies, arrive to similar results [37]. 
Among the many approaches pursued to evaluate the whole content of matter in the 
universe, there is the measurement of total mass relative to the baryon density. In this 
case we need the value of baryon density which implies the analysis of intracluster gas, 
the X-ray emission and the Sunyaev-Zeldovich effect (SZE). These measurements pro- 
vide Qm = 0.3 ±0.1 [38] consistent to the estimates obtained by different approaches. 
Tests performed to evaluate Qm are very important since they establish an upper limit 
to this parameter and give fundamental information in order to obtain indirect estimates 
of other parameters as in the cases of CMBR and SNela surveys. 

The SZE effect and the thermal bremsstrahlung (X-ray brightness data) for galaxy clus- 
ters represent intriguing astrophysical tests with several cosmological implications. In 
fact, distances measurements using SZE and X-ray emission from intracluster medium, 
are based on the fact that such processes depend on different combinations of some param- 
eters of clusters (e.g. see [39] and references therein). We recall that the SZE is a result 
of the inverse Compton scattering between CMBR photons and hot electrons of intra- 
cluster gas. The photon number is preserved, but photons gain energy and so a decrease 
of temperature is generated in the Rayleigh- Jeans region of black-body spectrum; on the 
other hand an increment is obtained in the Wien region. As it is well known, the Hubble 
constant H$ and the density parameters can be constrained by means of the angular 



d A (0, a)d A (a, a*) 
d A (0, a*) 



(13) 
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distances; it is possible to use the SZE and thermal bremsstrahlung to estimate such dis- 
tances. For a flat model with = 0.3, — 0.7, we have H = 63 ± 3kmsec~ l Mpc' 1 
while for an open universe with Qm = 0.3, it has been obtained Hq = 60 ± 3 [40]. 
There are many other approaches able to provide indications on cosmological parameters, 
however in this section we have just given a brief review of current status of observational 
results, avoiding of debating about errors, a relevant issue for which we remand to the 
literature. 

Besides, some of the above observational approaches and methods will be furtherly dis- 
cussed below in connection to the test of our models with observations. 



3 The Quintessence Approach. Are scalar fields strictly nec- 
essary to get acceleration? 

The above observational results lead to the straightforward conclusion that cosmic dy- 
namics cannot be explained in the traditional framework of standard model [9] so that 
further ingredients have to be introduced into the game. 

Several evidences suggest that, besides the four basic elements of cosmic matter-energy 
content, namely: baryons, leptons, photons and cold dark matter, we need a fifth element 
in order to explain apparent acceleration on extremely large scales. As we said above, 
under the standard of quintessence, we can enrol every ingredient capable of giving rise 
to such an acceleration. In other words, the old four element cosmology has to be sub- 
stituted by a new five elements cosmology. 

In this section, we will outline the new standard lore which was born in the last five or 
six years. The key element of such a new scheme is the fact that a scalar field can give 
rise to both the accelerated behaviour of cosmic fluid and the bulk of unclustered dark 
energy. 

Let us start from the cosmological Einstein-Friedmann equations which can be deduced 
from gravitational field equations when a FRW metric of the form (3) is imposed. We 
have: 

- = ~(p + 3p), (14) 
a 6 

d\ 2 k 1 , . 

aj a 1 3 

P + 3(£)(p + p) = 0, (16) 

where Eq.(14) is the Friedmann equation for the acceleration of scale factor a(t), Eq.(15) 
is the energy constraint and Eq.(16) is the continuity equation deduced from Bianchi 
contracted identities. We are using physical (Planck) units where 8irG — % — k B — 1 
unless otherwise stated. The source of these equations is a perfect fluid of standard 
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matter where p is the matter-energy density and p is the pressure of a generic fluid. 
A further equation is 

P = 7P, (17) 

which is the state equation which relates pressure and energy-density. In standard cos- 
mology 

0<7<1, (18) 
which is the so called Zeldovich interval where all the forms of ordinary matter fluid are 

enclosed (7 = gives "dust", 7 = - is "radiation", 7 = 1 is "stiff matter"). 

3 

c 2 

In standard units, it is 7 = -f where c s is the sound speed of the given fluid which cannot 

c 2 

exceed light speed. Introducing the relation (17), with the constraint (18), into Eq.(14), 
gives: 

p + 3p = p(l + 3 7 ) >0, (19) 

from which a < and then cosmological system cannot accelerate in any way. To get 
acceleration, the constraint (18) has to be relaxed so that a fluid of some exotic kind 
has to be taken into account. In other words, if the universe is found to be accelerating, 
there must be a matter-energy component pq with negative pressure pq such that 

Ptot + Sp tot < =s> p tot < - ^ < , (20) 

where p to t is the total matter-energy of the system. Obviously, 7 loses its physical inter- 
pretation of the squared ratio between speeds of sound and light. Since Pq + Pq > for 
any physically plausible negative pressure component and thanks to the positive energy 
conditions, pq > 0, we must have, at least, 

PQ> ?f^a>0, (21) 

where we are assuming all other components with non-negative pressure. In order to 
match all the above conditions, we use, as a paradigm, the evolution of a scalar field Q 
(quintessence scalar field) slowly rolling down its potential V(Q). In FRW space-time, 
we can define 

Pq = \(f + V(Q) , (22) 

which is the energy density of Q, 

PQ = \Q 2 ~ V(Q) , (23) 
which is the pressure. The evolution is given by the Klein-Gordon equation 

Q + 3HQ + V'(Q) = 0, (24) 
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while the Hubble parameter, in the case of a spatially flat space-time, is given by 



H 2 = = \(Pm + P Q ) , (25) 

where pu is the ordinary matter energy density. A suitable state equation is defined by 

= w = &-VIQ) 
Q PQ \Q 2 + V(Q)' 
and acceleration is obtained by the constraint 

-1 < 7 q < 0. (27) 

In this approach, quintessence is a time-varying component with negative pressure con- 
strained by Eq.(27). Vacuum energy density (i.e. cosmological constant) is quintessence 
in the limit 7q — > — 1. In order to give rise to structure formation quintessence could be 
also spatially inhomogeneous as some authors claim [11]. 

A crucial role is played by the potential V(Q) which has strict analogies with inflation- 
ary potentials. Also in the quintessence case, slow-roll is a condition in which the kinetic 
energy is less than the potential energy. Such a mechanism naturally produces a negative 
pressure. The key difference with inflation is that energy-scale for quintessence is much 
smaller and the associated time-scale is much longer compared to inflation. 
Several classes of potentials can give rise to interesting behaviours. For example, "run- 
away scalar fields" are promising candidates for quintessence. In this case, the potential 
V(Q), the slope V'(Q), the curvature V"(Q) and the ratios V'/V and V" /V have to 
converge to zero as Q — > oo (primes indicate derivative with respect to Q). 
The so-called tracking solution [11] are of particular interest in solving the coincidence 
problem. They are a sort of cosmological attractors which are obtained if the two main 
conditions 

V"V 2 

1q<1m r =~^T >! (28) 

holds. A specific kind of tracking solutions gives rise to the creeping quintessence which 
occurs for tracker potentials if initially pq ^> pu- Examples of potentials useful to get 
quintessence are 

V(Q) ~ Q- a , V(Q) ~ exp(l/Q) (29) 
but several other families are possible [12]. 

To conclude, we can say that a mechanism capable of producing accelerated expansion 
of cosmic flow can be obtained by introducing phenomenological scalar fields. As in the 
case of inflation, this is a paradigm which could be implemented also in other ways (for 
example, let us remember the Starobinsky inflationary scalaron [13]). 
The main criticism which could be risen to the above approach is that, up to now, no 
fundamental scalar field has been found acting as a quintessential dark energy field. For 
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example, there is no analogous of Higgs boson for quintessence. 

In the following, we propose other approaches with the aim to show that scalar fields are 
not strictly necessary to get quintessence. We will explore three possibilities: 

• The use of a better physically motivated equation of state (Van der Waals equation). 

• The extension of Einstein-Hilbert gravitational action including higher-order cur- 
vature invariant terms. 

• Considering torsion as a further source in cosmological equations. 



4 First hypothesis: Van der Waals Quintessence 
4.1 The model 

As previously said, the quintessential scheme can be achieved also by taking into account 
more general equations of state without adding any scalar field. The consideration we 
have done is the following. A perfect fluid equation of state is an approximation which 
is not always valid and which does not describe the phase transitions between successive 
thermodynamic states of cosmic fluids. In some epochs of cosmological evolution, for 
example at equivalence, two phases could be existed together. In these simple 
description by a perfect fluid equation of state could not be efficient. 

A straightforward generalization can be achieved by taking into account the Van der 
Waals equation of state which describes a two phase fluid. Also in this case, we have an 
approximation but it has interesting consequences on dynamics. 

Let us construct such a Van der Waals cosmology [41]. 

Instead of the equation of state p = jp, we take into account the Van der Waals 
equation [42] 

p-T^T p - ap '' (30) 

which gives the above one in the limit a, (5 — > 0. In standard units, the coefficients are 

a = 3p c vl = 3p c p~ 2 , (3 = \ c = (3p c ) _1 . (31) 

where p c , p c and v c are critical pressure, density and volume respectively. It is worth 
stressing the fact that we are using standard forms of matter as dust or radiation but the 
relation between the thermodynamic quantities p and p is more elaborate. The critical 
values are the indications that cosmic fluids change their phases at certain thermodynamic 
conditions. 

For the sake of simplicity, let us consider the case k = (FRW-flat cosmology) as 
strongly suggested by the CMBR data [7, 25, 44]. We have, from Eqs.(14), (15), (16): 



H 



-H 2 --(p + 3p), 



(32) 
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p = -3H(p + p), 



(33) 



H 2 



3 P - 



(34) 



The relation between p and p is given by Eq.(30), whereas Eq.(34) is a constraint, so 
that the effective variables of the system are {H, p] and the phase space is a plane. The 
singular points are given by the conditions 



H = p = 0, 

which, introduced into Eqs.(32) and (33), give 

p = -p. 

Immediately a de Sitter solution a(t) = ao exp Hot is found. Otherwise (H 
obtain the static solution a = a . 

The values of p which give such a situation are 



(35) 



(36) 
0), we 



Hi, 



Pl,2 

3 ' 



Pl,2 



(P + a)T yJ{P + «) 2 - Aa(3(l + 7) 



(37) 



2(1 + 7) 

so that we have two singular points at finite. They depends on the set of parameters 
{a,(3, , ~f} which, as we saw above, are connected to the critical thermodynamical pa- 
rameters {p c , v c , p c } of the fluid described by the Van der Waals equation. However, we 
discard the solutions 

#i, 2 = -^f, (38) 

which have no physical meaning in the today observed universe. 

In order to investigate the nature of such singular points, we have to "linearize" the 
dynamical system and study the local Lyapunov stability [43]. The system reduces to 



where 



A = 

B = 

C = 

D = 



H = AH + Bp, 
p = CH + Dp, 



Hl,2 , 

1 

+ 



7 



2 - «Pl,2 



6 2(l-/3 Pl , 2 ) 



-3|Pi,2 + - 1P l' 2 -2a Pl , 2 ) , 



-3#i, 



1 + 



1 - PPl,2 

7 



2ap 



1,2 



(39) 
(40) 



(41) 
(42) 

(43) 

(44) 
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where the constraints (30) and (34) have been taken into account. Since #1,2, Pi,2, ct, P, 
7 are real number, we have no imaginary component. The singular points can be nodes 
or saddle points. Their stability depends on the sign of A, B, C, D. If A, B,C,D < 0, 
the singular point is a stable node which attract the trajectories of phase space. If 
A, B,C,D > 0, the point is an unstable node. If the coefficients A, B, C, D are both 
positive and negative the singular point is a saddle. In the case of a stable node, the 
prescriptions of No Hair Conjecture are matched so that p\^ can be seen as effective 
cosmological constants. It is interesting to note that, in this case, such a cosmo logical 
constants are not imposed by hands but they are recovered from the critical values of the 
Van der Waals fluid, i.e. p c , v c , p c . 

The singular points at infinite can be studied by taking into account the equation 

dH = ^ 2 + |(p+ i ^-3V) 
dp 3tf( P+T ^-V) ' 

which is constructed by the dynamical system (32-34) and the Van der Waals relation 
(30). By imposing the asymptotic behaviour H ~ p m , we get H ~ p 1 / 2 as it has to 
be. However, prescriptions of No Hair Conjecture are recovered also in this case. It is 
interesting to observe that we get two values of cosmological "constant" at finite and at 
infinite which, in general, are different. 

The condition to obtain quintessence, i.e. an accelerated behaviour, is 

p + 3p < (46) 
which, in our specific case, corresponds to 




(47) 

In terms of the parameters of Van der Waals fluid, the conditions 

1 (3 + 3a ((3 + 3a) 2 n . 

have to hold at the same time in order to get a positive matter density. 

The observational constraints on are achieved in terms of parameters {o;,/3,7} 
being 



tt A = n A (a, [3, 7) = \ vz , 7 0.7 (49) 



which gives 



-l<g («,/3, 7 )<0. (50) 
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As usual, 

Pent = ^^2xl0- 29 gem- 3 (51) 
is the cosmological critical density. 

Such a qualitative discussion can be made quantitative by constraining Van der Waals 
cosmology with observations. 

4.2 Constraining Van der Waals quintessence by observation 

The three parameters (a, /3, 7) or (p c ,p c , 7) which feature Van der Waals scheme are not 
independent on each other. Considering the situation at critical points, one gets : 

3 , , 

Pc = g7Pc , (52) 

so that the number of independent parameters is now reduced to two, which are (/o c , 7). 
Using Eq.(52), we may rewrite Eq.(30) as: 

37^ 9 2 

V = 3^ - 8 7X (53) 
having introduced the scaled variables x = p/p c and y = p/p c - 

In this case, since in our model there is only one fluid, its nowaday density has to be 
equal to the critical density, i.e. : 

3H 2 1 

p(z = 0; p c , 7) = p crit -> p c = —— , (54) 

8itG £0(7) 

(hereinafter quantities labelled with refer to the today values, i.e. at redshift z — 0). 
At this point an important fact has to be stressed. It is clear that the r.h.s. of Eq.(53) 
can assume positive, negative and null values so that an "effective" matter-energy density 
could be defined. In order to remove such a "disturbing concept", we will generalize the 

c ^ 2 



parameter 7 which in our approach is not simply given by 7 = y — J but it fixes the 

relation between p e ff and p e //, two effective quantities assigned by the Van der Waals 
equation (30). With this consideration in mind 7 can assume negative values. 
Using Eqs.(52) and (54), we can completely characterize the cosmic fluid by only one 
parameter in the equation of state, in other words 7 become the only independent pa- 
rameter needed to describe the Van der Waals fluid. To constrain the value of 7, we can 
use some mathematical and physical considerations. Let us consider Eq.(16). We can 
rewrite it as a differential relation between the scaled density x and the scale factor a 
taking into account the Van der Waals relation: 

dx „ . I da „ da dx ,„ x 

_ = -3 (x + y) - ► -3 — = . (55) 

dt a dt a x + y 

It can be immediately integrated to give : 
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with : 



-3 (I+7) In (—) = In Pi (a;) + P 2 (x) arctan P 3 (x) 
\a J 



(56) 



Pi(x) 



X 



(57) 




^/9jx 2 - (277 + 8)x + 24(1 + 7) ' 
, II7-8 



(58) 




V1357 2 + 432 7 - 64 ' 
I87X - 277 - 8 



(59) 



V1357 2 + 432 7 - 64 ' 



Since the quantity on the left hand side of Eq.(56) is real, so must be the right hand 
side. We can use this simple mathematical condition to determine the range of significant 
values for 7. To this aim, it is useful to divide the real axis in three different regions 
bounded by the roots of the second order polynomial which is the function under the 
square root in P2(x) and Pz{x). Excluding the value 7 = — 1 in order to avoid the 
divergence of a due to the term (1 + 7) in Eq.(56), we have : 



We will examine the three different regions in Eq.(60) separately to see whether we can 
further constraint the range for 7 on the basis of physical considerations. 

4.2.1 Models with 7 e [ -3.34186, 0.141859 ] 

Let us start our analysis considering the second region in Eq.(60). For values of 7 in 
this range, the term under square root in Pi{x) and Pz{x) gets negative values, but the 
right hand side is still a real quantity. To show this, let us remember the two following 
mathematical relations \/b = iVb (b real and negative), arctanifr = i arctanhfr (b real), 
Eq.(56) may be rewritten as: 




00, -3.34186 [ U [ -3.34186, 0.141859 ] 



U ] 0.141859, 00 [ . 



(60) 



3(1+7) In 




) 



In Pi (x) + P^a^arctanhP^x) 



(61) 



with : 




II7-8 



(62) 




V-1357 2 -432 7 + 64 ' 
-l&yx + 277 + 8 



(63) 



V-1357 2 - 432 7 + 64 ' 



15 



xHz , g=-2 . OL 



• 1 


8 
6 












1 


4' 












1 


2 












-1 -0.5 

. 8 


. 5> 


^ 1 


1 . 5 


2 


2 . 5 



Figure 1: The scaled density x = p/p c as function of the redshift z for the model with 7 = —2.0. 
Note that x is confined between (x%, x^) (with 27(7 = —2.0) ~ 1.823 and ^2(7 = —2.0) ~ 0.731) 
as discussed in the text. 

Eq.(61) is correctly defined for 7, in the range we are considering, provided that the 
logarithmic term is real. To this aim, x should be in the range ] Xi, 22 [, being these 
latter the two roots of the square root term entering P\(x), given as: 

8 + 27 7 ±v/64 - 432 7 - 1357^ 
x 12 = . (64) 

7 

On the other hand, since x is positive by its definition, we must impose the constraint 
X\ t 2 > which means that we have to exclude all the values of 7 which lead to negative 
values of X\ or 2 2 - This constraint allows us to narrow the second region in Eq.(60) which 
is now reduced to : 

7 E [ -3.34186, -1 [ . (65) 

It is interesting to invert numerically Eq.(61) to get x(z; 7), being the red-shift z related 
to the scale factor as a /a = 1 + z. We plot the result for 7 = —2.0 in Fig. 1. This plot 
shows some interesting features. The scaled energy density x(z, 7) is bounded between 
two finite values as a consequence of the condition imposed so that the logarithmic term 
of Eq.(61) is correctly defined. It is interesting to observe that the model predicts that 
the energy density is an increasing function of time (and hence a decreasing function 
of red-shift as in the super-inflationary models). The two extreme values are (07, X2) 
that may be computed using Eq.(64) once a value of 7 in the range defined by Eq.(65) 
has been chosen. Furthermore, being x constrained in the range ]xi,X2[, we also get 
that the model may describe the dynamics of the universe over a limited red-shift range 
(z m in, z max ). This is not a serious shortcoming since the Van der Waals equation of state 
we are using is only an approximation (even if more realistic than the perfect fluid one) of 
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Figure 2: The scaled density x = p/p c as function of the redshift z for the model with 7 = 0.15. 
Note that the apparent divergence on the left side is not a physical one, but only an artifact of 
the scale used in the plot. 



the true unknown equation of state. It is thus not surprising that it may be applied only 
over a limited period of the universe evolution. It is quite easy to show (both analytically 
and numerically) that the lower limit of red-shift range is z m i n = — 1 so that the model 
may be used also to describe the near future evolution of universe. On the other hand, 
the upper limit is a function of 7, in particular, z max is always larger than z ~ 12 so 
that we can safely use the model to describe dynamics of the universe over the red-shift 
range (z < 1) probed by SNela Hubble diagram. We will see, in next subsection, that 
the situation will be different in the other red-shift ranges defined by Eq.(60). 

4.2.2 Models with 7 £ [ -3.34186, 0.141859 ] 

Let us turn now to the other two ranges defined in Eq.(60). It is quite easy to show that 
the term under square root, entering Pi(x) in Eq.(56), has the same sign of 7, if this 
parameter takes values in these regions. In order to have a positive function, we have to 
exclude the region ] — 00, —3.34186 [ which we do not consider anymore in the following 
so that we can now concentrate on the range ] 0.141859, 00 [. 

We can further narrow it on the basis of physical considerations. To this aim, let us 
invert numerically Eq.(56) to get x(z; 7). As an example, in Fig. 2 the results for 7 = 0.15 
are sketched. This plot shows two interesting features. First, there is a limiting red-shift 
zu m such that, for z > zu m , the density is not defined. This fact simply states that our 
model could describe the dynamical evolution of universe only for red-shifts lower than 
z Um . This fact could be explained qualitatively observing that x{zu m \ 7) = 3 for every 7. 
For p = 3p c , we get p — > 00. It is quite obvious that the universe cannot evolve until the 
pressure becomes finite (both positive or negative). Then, the universe starts evolving, 
but the density may follow two different tracks, decreasing or increasing with red-shift 
(that is increasing or decreasing with time). We will come back to this topic later. 
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We can observe that zu m is implicitly defined as function of 7 by the relation x(zu m ; 7) = 
3 and increases with 7 faster and faster as 7 approaches the lowest value in the range we 
are now considering. 

In the following section, we will investigate whether our model is able to reproduce the 
Hubble diagram of SNela. To this aim, the minimal request is that the model describes 
the dynamical evolution of matter density at least over the red-shift range probed by 
SNela, i.e. until z ~ 0.83. However, future observations (e.g. with the up to come SNAP 
satellite [47]) will be able to extend the upper limit of the range probed by SNela until 
z ~ 1.7. A supernova at this red-shift has been yet observed even if the estimated z is very 
uncertain and its apparent magnitude is likely to be seriously altered by gravitational 
lensing effects. 

Choosing a value of 7 very near to the lowest one will lead to zu m > 1000 so that our 
model could describe the universe from the decoupling age up to now. However, the Van 
der Waals equation of state is an approximation so that it should be not surprising if the 
model works only on a limited range for z. Thus, we impose an upper limit on the value 
of 7 by requiring that zu m {^i max ) = 1.7. This leads us to the following range for 7 : 



Before testing the models with 7 ranging in the interval defined by Eqs.(65) and (66), 
we come back again to Fig. 2. This plot shows that, for 7 e] 0.141859, 0.152370 [, the 
universe starts evolving when x 7^ 3, but there are two possible evolutionary tracks. 
From a mathematical point of view, this is a consequence of the intrinsic non linearity 
of Eq.(56). However, the two curves are clearly physically distinct since the density 
decreases or increases with the red-shift according to which track is followed. There is 
also another very important difference. Let us consider Eq.(14). If we ask for a today 
accelerating universe, we have to impose the following constraint on the present day value 
of the scaled density and pressure : 



For the values of 7 in the range defined by Eq.(66), the second degree polynomial in 
this equation is always positive definite, so that we get an accelerating universe if : 



Fig. 2 shows that the evolutionary track to get a today positive value of a is the upper 
one, so that we arrive to the conclusion that the accelerated models are those with an 
increasing density with time. Anyway, the result (from SNela) that the universe is now 
in an accelerated phase is model dependent so that, in the following, we will consider 
both models and test them against the observations to discriminate among the two 
evolutionary tracks. 



7 G ] 0.141859, 0.152370 [ . 



(66) 



xo + 3y < 




3 - x < 



x > 3 . 



18 



4.3 The Hubble diagram of SNela in the Van der Waals scheme 

It is well known that the use of astrophysical standard candles provides a fundamental 
tool to test different cosmological models. SNela are the best candidates to this aim since 
they can be accurately calibrated and can be detected at enough high red-shift. This fact 
allows to measure the SNela relation between magnitude and red-shift (i.e. the Hubble 
diagram) at high enough distances to discriminate among cosmological models. To this 
aim, one can match a given model with the observed Hubble diagram, conveniently 
expressed through Eq.(9), which we write as: 

/i(z) = 5\ogd L (z) + 25 , (67) 

being ji as before the distance modulus and D L (z) the luminosity distance calibrated 
with the Hubble constant. For our cosmological models, it is: 

d L (z) = ±- (1 +*) f [ P(z') Y 1 ' 2 dz> . (68) 

The distance modulus, as said, can be obtained from observations of SNela. The apparent 
magnitude m at the peak is measured, while the absolute one M may be deduced from 
template fitting or using the Multi Color Light Curve Shape method [4, 5]. The distance 
modulus is then /i = m — M. Finally the redshift z of the SNela is deduced from the host 
galaxy spectrum or (with a larger uncertainty) from the supernova spectrum directly. 
Our model can be fully characterized by two parameters : the today Hubble constant Ho 
and the value of 7. We find their best fit values minimizing the x 2 defined as 



rtJ/.. 7) -g ^fr *?■•'>-<•'' (69, 

i=l °M0>* ' ^mz,! 

where the sum is over the data points. In Eq.(69), a Mo is the estimated error on the dis- 
tance modulus, while a mz is the dispersion in the distance modulus due to the uncertainty 
a z on the measured red-shift. We have : 

Following Ref. [46], we assume a z = 200 kms -1 adding in quadrature 2500 km s _1 for 
those SNela whose red-shift is obtained from broad lines in their spectrum. Note that 
a mz depends on the parameters we wish to determine so that we have to use an iterative 
procedure to minimize the x 2 ■ 

The the SCP [4] and HZT team [5] have detected a quite large sample of high refshift 
(~ 0.18 0.83) SNela, while the Calan-Tololo survey [45] has investigated the nearby 
sources. Using the data in Refs. [4, 5], we have compiled a combined sample of 85 SNela 
as described in Ref. [46]. We exclude 6 likely outlier SNela as discussed in Ref. [4], thus 
ending with 79 SNela. 

Let us discuss first the results of the fit for the models with 7 e [ —3.34186, — 1 [. 
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Figure 3: 1, 2 and 3-cr confidence regions in the (Hq,j) plane for the models with 7 G 
] - 3.34186, -1 ]. 

As we can see in Fig.3 the models fit the data quite well. This is an encouraging 
result which mean that the proposed equation of state leads to a magnitude- red- shift 
relation which is in agreement with the observed SNela Hubble diagram. This result is 
further strengthened by the estimated value of the Hubble constant. The best fit value 
is H = 71 km s' 1 Mpc -1 , while the 68% confidence range turns out to be: 

H = (60 -=- 80) km s" 1 Mpc" 1 , 

which is in very good agreement with the most recent estimates in literature. For instance, 
the HST Key Project [49] has calibrated different local distance estimators and has found 
H = 72 ± 8 km s" 1 Mpc" 1 . 

It has to be stressed that the Hubble diagram of SNela (Fig.3) gives no constraint on the 
7 parameter. In fact, the best fit value for this range turns out to be 7 = —1.03, but the 
X 2 changes less than 1% for 7 in the range defined by Eq.(65). This is not an unexpected 
result since the study of density gives a density p(z, 7) almost completely independent of 
7 in the red-shift range probed by the available data [48]. 

However, the dependence on 7 becomes significative for higher red-shifts so that the 
prospects for using SNela to measure the value of 7 are quite good as soon as SNela at 
z ~ 1.0 1.7 (a red-shift range which will be explored by the dedicated SNAP satellite) 
will become available [47]. 

The results of the fit for models with 7 in the range given by Eq.(66) give also in this 
case that these models are able to reproduce the observed SNela Hubble diagram. We 
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Figure 4: 1, 2 and 3-cr confidence regions in the (Hq,j) plane for the models with 7 £ 
] 0.141859, 0.151237 ] having the scaled energy density increasing with the 5red-shift 

obtain H = 70 km s _1 Mpc -1 as best fit value for the Hubble constant, while the 68% 
confidence range is : 

H = (64 + 78) km s" 1 Mpc" 1 , 

which is still in very good agreement with the most recent estimates (see Fig. 4). However, 
the Hubble diagram of SNela gives no constraint on 7 parameter even for models with 7 
in the range defined by Eq.(66). The best fit value turns out to be 7 = 0.142, but as in 
the previous case it is not possible to give a reliable estimate of the uncertainty on this 
parameter. This is again due to the almost complete degeneracy of different models for 
z < 1 . 

On the other hand this degeneracy seems to be again eliminated using SNela data in 
the red-shift range z ~ 1.0 -j- 1.7. 

Finally, we have repeated the same analysis for models with density decreasing with 
time. Even if the fit is still possible, the estimated H is too low, the best fit value being 
40 km s -1 Mpc -1 . This fact has led us to do not consider anymore this class of models 
thus ending with the surprising conclusion that the models which fit the data are those 
with the energy density increasing with time. 

4.4 The age of the Universe in the Van der Waals scheme 

The SNela-Hubble diagram fit has shown that several Van der Waals models are able to 
reproduce the available data giving a value of H in very good agreement with the most 
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recent estimates. Furthermore, the matching with the data has allowed to reject models 
with p decreasing with time so that we do not consider them anymore in the following 
analysis. However, we are still not able to constrain the value of 7. To this aim, we 
can use the available estimates of the age of the universe. As a first step, we have to 
determine how the scale factor a depends on the time t. From Eq.(16) and the choice 
(53), we get : 



t-t lim = ^P- T [a'V^R)]- 1 da' , (71) 

-no Jaiim 

where au m is the value of a for the upper limit of red-shift range where the model may 
describe the evolution of the universe and t\ im is the corresponding time. We have also 
imposed a — 1 for the sake of shortness. This relation cannot be analytically inverted, 
but it is quite easy to invert numerically. Imposing a = 1 in Eq.(71), we have an estimate 
of to — tum, but we do not know the value of tu m , i.e. the age of the universe at a red-shift 
when the model starts to describe the dynamical evolution. On the other hand, we can 
evaluate the quantity : 



to-tref = ^r^[ 1 [a'v^]- 1 da', (72) 

tl Ja re f 

being t re f the age of the universe at a given red-shift z re f with a re f = a{z re f). This 
approach allows to constrain the values of 7, even if we do not know tu m , provided that 
we have an independent estimate of t re f. To this aim, we can use the age of an early-type 
high red-shift galaxy which can be estimated in a way which does not depend on the 
background cosmological model, but only on astrophysical assumptions on the galaxy 
formation scenario and stellar evolution. As a reference value, we use : 

z ref — 1-55 , t re f = ti.55 > 3.5 Gyr 

as first estimated in [50] for the radio galaxy 53W091 [51]. Different estimates of to are 
present in literature. Recent determinations of to from the age of globular clusters give 
t = 11.5 ± 1.5 Gyr [52, 53], while the analysis of the SNela Hubble diagram suggests t = 
14.9 (for H = 63km s~ l , Mpc~ r ) [4]. A multivariate analysis of combined BOOMERanG, 
DASI and MAXIMA data on the CMBR anisotropy spectrum gives estimates of t ranging 
from ~ 13 to ~ 16 Gyr according to the former ones used in the maximum likelihood 
estimation of the parameters [54]. We use t = 10 -7- 16 Gyr so we get : 

to - ti.ss = 6.5 -T- 12.5 Gyr (73) 

as a conservative limit on this quantity. 

For each value of 7 in the range defined by Eqs.(65) and (66) and for H in the 
corresponding 1 a confidence interval, we can determine the estimated range for to — ti.55. 
Imposing that this range is contained within the range given by Eq.(73), we find the 
constraint : 
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7 G [ -1.95, -1.40 ] (74) 

which we consider as our final estimate of the physically acceptable range for 7. Note 
that the range in Eq.(66) is completely cut out by this constraint since the predicted 
values of t — ti. 55 are much higher than the (quite conservative) limits we have imposed. 

4.5 Discussion of results in Van der Waals cosmology 

We have investigated the possibility that a cosmological model with a Van der Waals 
equation of state can give account of the accelerating expansion of the universe as sug- 
gested by the SNela data. 

As a first step, starting from mathematical and physical considerations, we have con- 
strained the effective parameter 7 of the model to take values in the range [ —3.14186, —1 [ 
U ] 0.141859, 0.15237 ]. In these regions, the matter energy density can be both increas- 
ing or decreasing with time as showed in Figs. 1 and 2 and the quantity x + 3y results 
always negative which is consistent with an accelerated phase of the universe evolution. 

As a following step, we have fitted these models with the data of SNela-Hubble dia- 
gram. The fitting procedure has lead us to conclude that both the intervals of 7 are 
physically significative, but the decreasing behaviour of energy density in the range 
]0. 141859, 0.15237] has to be discarded since the estimated Hubble constant is only 
40 km s" 1 Mpc -1 , in strong disagreement with results in literature. On the other hand, 
the best fit values of the parameters turn out to be : 

7 = -1.09 , H = 71 km s" 1 Mpc" 1 , 
for models with 7 G ] — 3.34186, —1 [ and: 

7 = 0.142 , H = 70 km s' 1 Mpc" 1 , 

for models with 7 G ] 0.14186, 0.15237 [ and the energy density increasing with time. 
The results for Hq are in perfect agreement with the estimates obtained with different 
methods such as the calibration of local distance estimators [49], which is an evidence of 
the validity of our models. Unfortunately, the very small variation of the energy density 
with 7 in the red-shift range probed by available data does not let us to discriminate in 
a useful way between different values of this parameter. However, higher red-shift SNela 
data (as the ones that will be collected by the next to come SNAP satellite) will be able 
to break the degeneracy among the different values of 7 since they will probe a red-shift 
region where the energy density is more sensitive to the effective parameter of the theory. 

As a possible test to break the degeneracy among the values of 7 with nowaday 
available data, we have estimated the quantity t Q — £1.55, being t the nowaday age of the 
universe and ti.55 = t(z = 1.55). For a fixed 7, we estimate a range for t — £1.55 changing 
Hq in the corresponding 68% confidence region and compare it with the conservative 
limit given in Eq.(73). This test allows us to exclude the range in Eq.(66) since the 
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predicted values of to — £1.55 are much higher than 12.5 Gyr so that we turn out with a 
final estimate for 7 as : 

7 E [ -1.95, -1.40 ] . 

The analysis presented suggests that a cosmological model characterized by a Van der 
Waals equation of state is a novel and interesting approach to explain the accelerating 
expansion of the universe. Indeed, for values of 7 in the range given by Eq.(74), the 
models are able to fit quite well the SNela Hubble diagram and predict an age of the 
universe consistent with recent estimates obtained with completely different methods. 
Another feature of the selected models is the behaviour of the energy density which is an 
increasing function of time (i.e. decreasing with z as shown in Fig. 1) bounded between 
two finite values. However, some models of super- inflation also predict such a behaviour 
so that it should be interesting to see whether a cosmological fluid with a Van der Waals 
equation of state could be accomodated within the framework of these theories. 

5 Second hypothesis: Curvature Quintessence 
5.1 The model 

There is no a priori reason to restrict the gravitational Lagrangian to a linear function of 
Ricci scalar R minimally coupled with matter [55]. Furthermore, we have to note that, 
recently, some authors have taken into serious consideration the idea that there are no 
"exact" laws of physics, but that Lagrangians of physical interactions are "stochastic" 
functions with the property that local gauge invariances (i.e. conservation laws) are 
well approximated in the low energy limit with the properties that fundamental physical 
constants can vary [56]. This scheme was adopted in order to treat the quantization on 
curved spacetimes. The result was that either interactions among quantum fields and 
background geometry or the gravitational self-interactions yield corrective terms in the 
Einstein-Hilbert Lagrangian [57]. Futhermore, it has been realized that such corrective 
terms are inescapable if we want to obtain the effective action of quantum gravity on scales 
closed to the Planck length [58]. They are higher-order terms in curvature invariants as 
R 2 , R^R^v, R^^R^ap, ROR y or RO k R y or nonminimally coupled terms between scalar 
fields and geometry as 4> 2 R. Terms of this kind arise also in the effective Lagrangian 
of strings and Kaluza-Klein theories when the mechanism of dimensional reduction is 
working [59]. 

Besides fundamental physics motivations, all these theories have acquired a huge 
interest in cosmology due to the fact that they "naturally" exhibit inflationary behaviours 
and that the related cosmological models seem very realistic [60, 61]. Furthermore, it is 
possible to show that, via conformal transformations, the higher-order and nonminimally 
coupled terms {Jordan frame) always corresponds to the Einstein gravity plus one or more 
than one minimally coupled scalar fields {Einstein frame) [62, 63, 64, 65, 66] so that these 
geometric contributions can always have a "matter-like" interpretation. 
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Here we want to investigate the possibility that quintessence could be achieved by extra- 
curvature contributions. We focus our attention on fourth-order theories of gravity [67]. 
A generic fourth-order theory in four dimensions can be described by the action 

A = J d 4 x^g~[f(R) + L M } , (75) 

where f(R) is a function of Ricci scalar R and Lm is the standard matter Lagrangian 
density. We are using physical units 8tiGn — c — K — 1. The field equations are 

f(R)Rcfi - \f{R)g aP = f'(Ry aP (g«»gp» - g*^) + , (76) 

which can be recast in the more expressive form 

G a p = R a /3 - -g a pR = T£p V + , (77) 

where 

?T = j^j {\g«P \f(R) - Rf'(R)\ + f'(Ry> a/3 (ga^ - g*^)} (78) 

and 



rpM _ _J rpM , -n i 

13 ~ f'(R) 13 ' 

is the stress-energy tensor of matter where we have taken into account the nonmnimal 
coupling to geometry. The prime indicates the derivative with respect to R. If f(R) = 
R + 2A, the standard second-order gravity is recovered. 

It is possible to reduce the action to a point-like, FRW one. We have to write 

A CU rv = J dtC(a, a; R, R) , (80) 

where dot indicates derivative with respect to the cosmic time. The scale factor a and 
the Ricci scalar R are the canonical variables. This position could seem arbitrary since 
R depends on a, a, a, but it is generally used in canonical quantization [68, 69, 70]. The 
definition of R in terms of a, a, a introduces a constraint which eliminates second and 
higher order derivatives in action (80), and gives a system of second order differential 
equations in {a, R}. Action (80) can be written as 



A curv = 2tt 2 J dt | a 3 f(R) - A 



„ (da k 
R + 61- + — + — 



(81) 



where the Lagrange multiplier A is derived by varying with respect to R. It is 

A = a 3 f'(R) . (82) 
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The point-like Lagrangian is then 

C = C curv + C M = a 3 [f(R) - Rf'(R)} + Qaa 2 f(R) + 

+Qa 2 aRf"(R) - Qkaf'(R) + a 3 p M , (83) 

where we have taken into account also the fluid matter contribution which is, essentially, 
a pressure term [71] . 

The Euler-Lagrange equations are 

a\ / a\ 2 k , . 

+ (-) +^ = -Ptot, (84) 



and 



f"(R) 



ay \ay a 



„ (a a 2,2 k s 
R + 61- + — +^ 
\a a a z , 



. (85) 



The dynamical system is completed by the energy condition 



d\ 2 k 1 
aj a z 3 



Combining Eq.(84) and Eq.(86), we obtain 

=-\^t + ^Ptot] , (87) 

where it is clear that the accelerated or decelerated behaviour depends on the r.h.s. How 
it is evident equations (84)-(87) generalize in a very simply way the system (14)-(16) for 
standard cosmological model. 
Now it holds: 

Ptot — Pcurv + Pm Ptot — + Pm , (88) 

where we have distinguished the curvature and matter contributions. 

In fact by curvature-stress-energy tensor, we can define a curvature pressure 

P» = j^j {2 (£) Rf"(R) + Rf"(R) + R 2 f"{R) - \ [f(R) - Rf'(R)]} , (89) 

and a curvature density 

Pour, = j^j {\ lf(R) - Rf'(R)} - 3 (£) Rf"(R)} . (90) 

From Eq.(87), the accelerated behaviour is achieved if 

Ptot + Zptot < , (91) 
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which means 

Pcurv > T^Ptot , (92) 

assuming that all matter components have non-negative pressure. 

In other words, conditions to obtain acceleration depends on the relation 



Pcurv + S Pcurv = [R 2 f"\R) + (£) Rf\R) + Rf"(R) - I [f(R) - Rf'(R)}) , 

(93) 

which has to be compared with matter contribution. However, it has to be 

P( - curv ^ = ~ -Ki < (94) 

Icurv j - 1 - _ Icurv ^ u • W / 

Pcurv 

The form of f(R) is the main ingredient to obtain this curvature quintessence. As simple 
choice in order to fit the above prescriptions, we ask for solutions of the form 

f{R) = UR\ a(t) = aoQ-y . (95) 

The interesting cases are for n ^ 1 (being n = 1 Einstein gravity) and a > 1 (accelerated 
behaviour). Inserting Eqs.(95) into the above dynamical system, we obtain the exact 
solutions 

a = 2; n= -1,3/2; k = . (96) 
In both cases, the deceleration parameter is 

Qo = ~, (97) 

in perfect agreement with the observational results. 

The case n = 3/2 deserves further discussion. It is interesting in conformal transfor- 
mations from Jordan frame to Einstein frame [72, 73] since it is possible to give explicit 
form of scalar field potential. In fact, if 



9ap = f{R)g a p, y» = ^|ln/'(i2), 
we have the conformal equivalence of the Lagrangians 

£=V=gfoR 3/2 ^ 1= xHi 



(98) 



Y + ^y?VV - ^oexp I J^cp 



(99) 



in our physical units. This is the so-called Liouville field theory and it is one of the few 
cases where a fourth-order Lagrangian can be expressed, in the Einstein frame, in terms 
of elementary functions under a conformal transformation. Furthermore this is a case in 
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which quintessence can be achieved also by an exponential potential. 
It is possible to obtain the general solution [74] 

a(t) = a [at 4 + c 3 t 3 + c 2 t 2 + at + c ] 1/2 . (100) 



The constants a are combinations of the initial conditions. Their particular values deter- 
mine the type of cosmological evolution. For example, a 7^ gives a power law inflation 
while, if the regime is dominated by the linear term in a, we get a radiation-dominated 
stage. 

These cases belong to an entire family of exact solutions for the quintessence curvature 
model. In fact solving Eqs. (87) and (88) in the limit in which matter can be neglected 
(this hypothesis can be considered well posed in a toy model) and assuming for f(R) and 
a(t) behaviours like (95), we obtain the algebraic system for the parameters n and a 

' a[a{n - 2) + 2n 2 - 3n + 1] = 

< ' (101) 
k a[n 2 + a(n - 2 - n - 1)] = n(n - l)(2n - 1) 



from which we have the solutions: 

a = — > n = 0, 1/2, 1 

_2n 2 + 3n-l (102) 

a = , n . 

n-2 



The cases for a = give static models not interesting to discuss. The cases for a and 

n generic deserve more attention since a class of cosmological quintessential behaviours 

can be found. In Fig. 5 there is a sketch of the situation. 

mi r 1 1 c 1 • 1 —2n 2 + 3n — 1 . 

ihe state equation for the class of solutions n and a = gives 

n — 2 

Gn 2 -In -I 

7 - = - 6n 2 -9n + 3 " (1 ° 3) 



Accelerated behaviours for increasing scale factors are possible only for 7 cm ,„ < how 
requested for a cosmological fluid with negative pressure. The whole approach seems 
intriguing in relation to the possibility to get an accelerated phase of universe expansion 
as an effect of a higher-order curvature Lagrangian. 

At this point, to obtain constraints on n and to match with the observations we take 
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20 - 




-40 - 

Figure 5: The plot shows the behaviour of a and against n. It is evident a region in which 
the power law of the scale factor is more than one and correspondently the parameter 7( cnri) ) 
of state equation is negative . 



into account data, in particular we compare this model in relation to SNela results as 
yet done for the Van der Waals case in ([48] and [75]). To complete our analysis we have 
to evaluate also the age prediction of curvature quintessence model as a function of the 
parameter n. 

In this sense, some results are present in literature where the form of f(R) is selected by 
the CMBR constraint [76]. 

5.2 Matching Curvature Quintessence with observations 
5.2.1 SNela matching 

In order to verify if such an approach has real perspective to be physically acceptable, we 
match the theory with observational data. In this way, we can constraint the parameters. 
To this aim we consider the supernovae observations, in particular we use the data in 
[4, 5] as in the Van der Waals case. 

Starting from these data, it is possible to perform a comparison between the theoretically 
predicted expression of distance modulus and its observational value. In this case the 
luminosity distance can be expressed by the general expression: 




where E(() = — and c is the light speed having reintroduced the standard units. 
H 
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Starting from Eq. (95), the Hubble parameter, in terms of redshift, is given by the 
relation 



H(z) = H (l + z) 1/a (105) 

where a is given in (102) as a function of n. 
The luminosity distance is now: 

d L (z,H ,n) = ^-(l + z) [\l + ()-y«d( (106) 
-no Jo 

and after the integration, we obtain 

C / OL \ r a i 

d L (z, H , n) = — ( (1 + z) (1 + z)— - 1 . (107) 

H \a — 1/ L J 

This expression is not defined for a — 0, 1, so we have to check the luminosity distance 

taking into account the existence of such singularities. As a consequence the fit can be 

1 1 11 

performed in five intervals of n, that is: n < -(1 — v^3) , -(1 — V^3) < n < - , - < n < 

Z Z Z Z 

1 , K n< \(1 + V3) , n > \{l + VS). 
z z 

Hubble parameter, as a function of n, shows the same trend of a (Fig.5). We find that for 
n negative, lower than -100, it is strictly increasing while for n positive, greater than 100, 
it is strictly decreasing. In relation to this feature, we have tested n for values ranging 
in these limits because, as we shall see below, outside of this range the value of the age 
of universe becomes manifestly not physically significant. The results can be showed in 
the Table 1. From Table 1, in relation to the best fit values of H and x 2 > we can only 

exclude the range - < n < 1 as physically not interesting. In the other cases, the results 

give interesting best fit values both for the Hubble parameter and the x 2 , indicating that 
such a model could represent a significant theoretical background to explain SNela data. 



We have to note that the best fit is completely degenerate with respect to the n 
parameter. This peculiarity indicates that x 2 varies slightly with respect to n hindering 
the possibilities of constrain n. 
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Range 


H% est (km sec~ 1 Mpc~ 1 ) 






-100 < n < 1/2(1 - y/3) 


65 


-0.73 


1.003 


1/2(1 - V3) < n < 1/2 


63 


-0.36 


1.160 


1/2 < n< 1 


100 


0.78 


348.97 


1< n < 1/2(1 + v 7 ^) 


62 


1.36 


1.182 


1/2(1 + V3) <n<3 


65 


1.45 


1.003 


3 < n < 100 


70 


100 


1.418 



Table 1: Results obtained fitting the curvature quintessence model with Snla data. First 
column indicate the range of n studied, column two gives the relative best fit Ho, column three 
n best , column four the \ 2 index. 

5.2.2 The age of universe in Curvature Quintessence approach 

The age of the universe can be simply obtained from a theoretical point of view if one 
has the value of the Hubble parameter. Now, in curvature quintessence model from the 
definition of H, using the relation (95), we have: 

t = aH' 1 (108) 

which from a (102) becomes 

^ _ 2 „, + 3 „- 1 /r , (iog) 

so, the today value of the age of universe simply is t — a^H^ 1 

By a simple algorithm we can evaluate the age of universe taking into account both 
the intervals and the la range of variability of the Hubble parameter, deduced from the 
supernovae fit for each interval. We have considered as good predictions, age estimates 
included between 10 and 18 Gyr. 

First of all we have discarded the intervals of n which give a negative value of t. Obvi- 
ously from Eq.(108), negative values for the age of the universe are obtained for negative 
values of a, so we have to exclude the ranges 1/2 < n < 1 and n > 2 (Fig.5). The results 
are shown in Table 2. 



As we have said above this test gives interesting results. 

The last check for our model is to verify if the range of n significant for the data provides 
also an accelerated rate of expansion. This test can be easily performed starting by re- 
lations (4), (95) and (102). To have an accelerated behaviour, the scale factor function 
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Range 


AH(km sec 1 Mpc 1 ) 


An 


t(n best )(Gyr) 


-100 < n < 1/2(1 - V3) 


50-80 


-0.67 < n < -0.37 


23.4 


1/2(1 - v^) < n < 1/2 


57-69 


-0.37 < n < -0.07 


15.6 


1 < n < 1/2(1 + V3) 


56-70 


1.28 < n < 1.36 


15.3 


1/2(1 + v^3) < n < 2 


54-78 


1.37 < n < 1.43 


24.6 



Table 2: The results of the age test. In the first column is presented the tested range. Second 
column shows the la Hq range obtained by supernovae test, while in the third we give the n 
intervals, that is the values of n which allow to obtain ages of the universe comprised between 
10 and 18 Gyr. In the last column, the values of the age obtained for the best fit value of each 
interval are reported. 

a = a t a has to get values of a negative or positive greater than one. We obtain that in- 
tervals —0.67 < n < 0.37 and 1.37 < n < 1.43 provide a negative deceleration parameter 
with a > 1, that is cosmological models expanding with an accelerated rate. Conversely 
the other two intervals of Table 2 do not give interesting cosmological dynamics in this 
sense (q > 0, < a < 1). 

Finally we are able to state that cosmological models based on a relativistic lagrangian 
with a generic power of R (different from one) are physically coherent with observations. 
This results indicate again that quintessence can be achieved without scalar fields starting 
from effective fundamental theories. 

6 Third hypothesis: Torsion quintessence 
6.1 The model 

Taking into account torsion is a straightforward generalization to implement concepts 
as spin in General Relativity [77, 78]. However, it was soon evident that torsion, as 
considered e.g. in Einstein-Cartan-Sciama-Kibble (ECKS) theory, does not seems to 
give relevant effects in the observed astrophysical structures. Nevertheless it was found 
that for densities of the order of 10 47 g/cm 3 for electrons and 10 54 g/cm 3 for protons and 
neutrons, torsion could give observable consequences if all the spins of the particles are 
aligned. These huge densities can be reached only in the early universe so that cosmology 
is the only viable approach to test torsion effects [79]. However no relevant tests con- 
firming the presence of torsion have been found until now and it is still an open debate 
if the space-time is Riemannian or not. Considering the cosmological point of view and, 
in particular the primordial phase transitions and inflation [1, 2, 80], it seems very likely 
that, in some regions of the early universe, the presence of local magnetic fields could 
have aligned the spins of particles. At very high densities, this effect could influence the 
evolution of primordial perturbations remaining as an imprint in today observed large 
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scale structures. In other words, a main goal could be to select perturbation scales con- 
nected to the presence of torsion in early epochs which give today-observable cosmological 
effects [81]. 

From another point of view, the presence of torsion could give observable effects with- 
out taking into account clustered matter. Here, we want to investigate if the quintessen- 
tial scheme can be achieved by taking into account theories of gravity with torsion [82]. 

In the ECSK theory, the affine connection is non-symmetric in its lower indices and 
the antisymmetric part 

rf &c] = s bc a (no) 

is called torsion. Usually, we can divide such an antisymmetric part into three compo- 
nents: one of them is irreducible, while the other two can be set to zero [83, 85]. This 
assumption yields the simplest theory containing torsion. 
Furthermore, we can express by a 4- vector 

a a = e abcd Sbcd (m) 

the totally antisymmetric part of torsion. If one imposes to it the symmetries of a 
background which is homogeneous and isotropic, it follows that, in comoving coordinates, 
only the component <7° survives as a function depending only on cosmic time (see [83, 
85, 86]). 

For a perfect fluid, the Einstein-Friedmann cosmological equations in presence of 
torsion can be written as usual 

^ = 4(p + 3p), (112) 

and 

'a\ 2 k 1_ 

- +- = ttP- (113) 

.a J cr 3 

Energy density and pressure can be assumed in the forms [83, 85] 

~P = P + f, P = P~f, (114) 

where / is a function related to a , while p and p are the usual quantities of General 
Relativity. This choice can be pursued since we can define S a b c = S[ a b c ] = fit) where f(t) 
is a generic function of time which we consider as the source of torsion. For a detailed 
discussion of this point see [85]. 

As usual, we can define a stress-energy tensor of the form 

= (p + p)u a u b - pg ab , (115) 
which, by Eqs. (114), can be splitted as 

rptot rpM I rptorsion /l 1M 

1 ab — 1 ab + 1 ab ■ l ii0 J 
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Due to the contracted Bianchi identity, we have 

T^ ;b = 0; (117) 
from which, we can assume that (cfr. [84]) 

TM;b n rptorsion ;b n /1 1 o\ 

ab = U > 1 ab =U - 

In the FRW space-time, Eq.(117) becomes 

P + 3H(p + p) = (119) 

which is 

p + 3H(p + p) = -2ff. (120) 

From Eqs.(118), both sides of (120) vanish independently so that 

f(t) = f = constant. (121) 

In other worlds, a torsion field gives rise to a constant energy density, i.e. a torsion — A 
term. Taking into account standard matter which equation of state is defined into the 
above Zeldovich range, we obtain 



P = Po 



3(7+1) 

+ fo ■ (122) 



Oo 

. a _ 

Inserting this result into the cosmo logical equations, we get, in any case, a monotonic 
expansion being /q > 0, a 2 > 0. The condition to obtain the accelerated behaviour is 

p + 3p< 2/ 2 , (123) 

so that acceleration depends on the torsion density. 
In a dust-dominated universe, we have 

P = Po(^) 3 + / 2 , P = -fl (124) 
and then solving the cosmological equations we obtain the general solution: 

a(t)= (|^) 1/3[cOSh(/ot) - 1]1/3 - (125) 
Obviously, if fot — > 0, we have cosh(/ £) — 1 — (fot) 2 and then a ~ t 2 / 3 , as it has to 

be. 

This results tell us that a sort of dark energy (a torsion-A term) can be obtained 
without considering additional scalar fields in the dynamics but only assuming that the 
space-time is U4 instead of V4. A preminent role is played by the type of torsion we are 
going to consider [85]. However, such a density should be comparable to the observed 
limits of dark energy (i.e. ~ 0.65 -h 0.7) in order to give relevant effects. Another 
point is that such a torsion quintessence should match the issues of cosmic coincidence 
[11] as scalar field quintessence. This point strictly depends on the type of torsion since 
torsion can be or not related to the spin density [85]. In the first case, the spin of baryonic 
and non-baryonic matter would rule the dark energy (torsion) density. 
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6.2 Matching Torsion quintessence with observations 



6.2.1 The Supernovae SNe la method 

As in previous cases we compare the theoretical estimates of luminosity distance with 
the results obtained for type la supernovae. Our analysis is based again on the sample 
data [46] from SCP and HZT observations [4, 5] 

The luminosity distance in the model we are considering is completely equivalent to 
the one in a spatially flat universe with a non-zero cosmological constant. Thus <1l{z) is 
simply given as : 

d L (z) = (l + z) f dz'[n M (l + z'f + Qtor]' 1 / 2 . (126) 
Jo 

where f2 tor = 1 — plays the same role as the usual Q\. 

Our model can be, so, fully characterized by two parameters : the today Hubble 
constant H and the matter density Qm- We find their best fit values minimizing the x 2 
defined in this case as : 

X 2 (^o,M=E rr( f ^ ) ' /ifS12 (127) 



where the sum is over the data points [46]. 

We remember that the <t Mo and a rnz are the errors related to the distance modulus eval- 
uation, their estimate is done as in the previous cases [46]. 

The results of the fit are presented in Fig. 6 where we show the 1,2 and 3 c confidence 
regions in the (Qm, H ) plane. 

The best fit values (with la error) turn out to be : 

Vt M = 0.30 ±0.08 , H = 65.4 ± 1.2 km s' 1 Mpc' 1 . 

which allow to conclude that a torsion A-term could explain observation very well. On 
the other hand, we can estimate the torsion density contributions which result to be : 

fl = (5.6 ± 0.7) x 1(T 30 g cm' 3 , 

which is a good value if compared to the cosmological critical density. It is worthwhile 
to note that, in the case of SNela, the error on H does not take into account systematic 
uncertainties due to possible calibration errors. 



6.2.2 The Sunyaev-Zeldovich/X-ray method 

Besides the above results, we can discuss how the Hubble constant H and the torsion 
density parameter Q tor can be constrained also by the angular diameter distance Da 
as measured using the Sunyaev-Zeldovich effect (SZE) and the thermal bremsstrahlung 
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(X-ray brightness data) for galaxy clusters. 

We will limit our analysis to the so called thermal or static SZE, which is present in 
all the clusters, neglecting the kinematic effect, which is present in those clusters with a 
nonzero peculiar velocity with respect to the Hubble flow along the line of sight. Typically 
the thermal SZE is an order of magnitude larger than the kinematic one. The shift of 
temperature is: 

(128) 

where x = - — — is a dimensionless variable, T is the radiation temperature, and y is the 
k B T 

so called Compton parameter, defined as the optical depth r = <jt J n e dl times the energy 
gain per scattering: 

y = [ ^n e a T dl. (129) 
J m e c 2 

In Eq. (129), T e is the temperature of the electrons in the intracluster gas, m e is the 
electron mass, n e is the numerical density of electrons, and a T is the cross section of 
Thompson electron scattering. We have used the condition T e 3> T (T e is of the order 
10 7 K and T, which is the CMBR temperature is ~ 2.7 K). Considering the low frequency 



iXl , fx . 

— = y rrcoth I - ) - I 
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regime of the Rayleigh- Jeans approximation, we obtain 



T 



2y (130) 



The next step to quantify the SZE decrement is that we need to specify the models for 
the intracluster electron density and temperature distribution. The most commonly used 
model is the so called isothermal (3 model [87] . We have 



2\ ~¥ 



n e (r) = n e (r) = n eo |l + \-j J , (131) 

T e (r) = T eo , (132) 

being n eo and T eo , respectively the central electron number density and temperature of 
the intracluster electron gas, r e and f3 are fitting parameters connected with the model 
[88]. For the effect of cluster modelling see [89]. From Eq. (129) we have 

Jo V m e c z ) 

being 

E ^r( i+ ©y T *- (i34) 



The integral in Eq. (134) is overestimated since clusters have a finite radius. The effects 
of the finite extension of the cluters are analyzed in [89, 90]. 

A simple geometrical argument converts the integral in Eq. (134) in angular form, by 
introducing the angular diameter distance, so that 

/ /tf\ 2 \ 1/2-3/3/2 r(&-±) 

E = 6> c 1+ - ^4kr^ (135) 



02 



r(f) 



c 



In terms of the dimensionless angular diameter distances, cIa (such that Da = -rrdA) 

Ho 



we get 



30 lW / „ \ 2\ hi 1 " 3 ^ 



AT(0) 2 ( a T K B T ec n eo \ ^ [fj (9, , 

—=-h- { m eC J^Y^yy-yrJ I dA > (136) 

and, consequently, for the central temperature decrement, we get 



AT(0 = O)_ 2 fa T K B T ec n eo \ r(¥§) 



T H V m e c J v r (f) Ho 



c 



-^T^irdA. (137) 
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The factor — <Ia in Eq. (137) carrys the dependence on the thermal SZE on both the 

cosmological models (through H and the Dyer-Roeder distance (La) and the redshift 
(through cLa)- From Eq. (137), we can note that the central electron number density is 
proportional to the inverse of the angular diameter distance, when calculated through 
SZE measurements. This circumstance allows to determine the distance of cluster, and 
then the Hubble constant, by the measurements of its thermal SZE and its X-ray emission. 
This possibility is based on the different power laws, according to which the decrement 

A'T(f) — ft) 

of the temperature in SZE, — , and X-ray emissivity, Sx, scale with respect to 

To 

the electron density. In fact, as above pointed out, the electron density, when calculated 
from SZE data, scales as d~ A x ( nf ZE oc d^ 1 ), while the same one scales as d~ A 2 (nfo~ ray oc 
d A 2 ) when calculated from X-ray data. Actually, for the X-ray surface brightness, Sx, 
assuming for the temperature distribution of T e = T e0 , we get the following formula: 

Sx = % n2 eoj^y 3 0c^- Q dAlsx, (138) 

being 

h x = r (—) 2 di, 

Jo \n e oJ 

the X-ray structure integral, and ex the spectral emissivity of the gas (which, for T e > 
3xl0 7 , can be approximated by a typical value: ex = £\/T~ e , , with e ~ 3.0xl0~ 27 np erg 
cm" 3 s" 1 K^ 1 [88]) . The angular diameter distance can be deduced by eliminating the 
electron density from Eqs. (137) and (138), yielding: 



y 2 An(l + z 




(139) 

\m e u- j no u\Ojj — g, 2) 

where B(a, b) = — is the Beta function. 
T(a + b) 

It turns out that 

D A = ±4 A cc ^1 (140) 

-"0 ^>X0 1 eO U c 

where all these quantities are evaluated along the line of sight towards the cluster center 
(subscript 0), and 9 C is referred to a characteristic cluster scale along the line of sight. It 
is evident that the specific meaning of this scale depends on the density model adopted 
for the cluster. In our calculations we are using the so called (3 model. 

Eqs. (139) allows to compute the Hubble constant H , once the redshift of cluster is 
known and the other cosmological parameters are, in same way, constrained. Since the 
dimensionless Dyer-Roeder distance, dA, depends on fltor, ^m, comparing the estimated 
values with the theoretical formulas for Da, it is possible to obtain information about Qm 
ttton an d Hq. Recently, distances of 18 clusters with redshift ranging from z ~ 0.14 to z ~ 
0.78 have been determined from a likelihood joint analysis of SZE and X-ray observations 
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(see [91] and reference therein). Modeling the intracluster gas as a spherical isothermal 
/3-model allows to obtain constraints on the Hubble constant H in a standard A-FRW 
model. We perform a similar analysis using angular diameter distances measurements 
for a sample of 44 clusters, constituted by the 18 above quoted clusters and other 24 
already know data (see [39]). 

As indicated in [39, 91], the errors a are only of statistical nature. Taking into account 
our model with torsion (125), the theoretical expression for the angular diameter distances 
D A is 

D A (z) = —±—d L (z) = -L_ [* dz'[n M (l + z'f + n tor Y 1 ' 2 . (141) 
[L + z) 1 + z Jo 

We find the best fit values for f2 tor and H , minimizing the reduced x 2 : 



x\h , n M ) = £ \ D ^ theor ^fM)-D A tr ^ (142) 



The best fit values (at la) turn out to be : 

tt M = 0.30 ±0.3 , H = 68 ± 6 km s' 1 Mpc' 1 
in good agreement with the above fit derived from SNe la data. 



In conclusion, we have shown that the net effect of torsion is the introduction of an 
extra-term into fluid matter density and pressure which is capable of giving rise to an 
accelerated behaviour of cosmic fluid. Being such a term a constant, we can consider it a 
sort of torsion A-term. If the standard fluid matter is dust, we can exactly solve dynamics 
which is in agreement with the usual Friedmann model (to be precise Einstein-de Sitter) 
as soon as torsion contribution approaches to zero. 

The next step has been to compare the result with observations in order to see if 
such a torsion cosmology gives rise to a coherent picture. We have used SNe la data, 
Sunyaev-Zeldovich effect and X-ray emission from galaxy clusters. Using our model, we 
are capable to reproduce the best fit values of H and Qm which gives a cosmological 
model dominated by a cosmological A-term. In other words, it seems that introducing 
torsion (and then spins) in dynamics allows to explain in a natural way the presence 
of cosmological constant or a generic form of dark energy without the introduction of 
exotic scalar fields. Besides, observations allows to estimate torsion density which can 
be comparable to other forms of matter energy (~ 5. 5x 10~ 30 g cm -3 ). 

However, we have to say that we used only a particular form of torsion and the 
argument can be more general if extended to all the forms of torsion [85]. Furthermore, 
being in our case the torsion contribution a constant density, it is not possible to solve 
coincidence and fine tuning problems. To address these issues we need a form of torsion 
evolving with time. 
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7 Discussion, conclusions and perspectives 

In this review paper, we have discussed, from theoretical and observational points of view, 
some approaches by which it is possible to obtain the observed accelerated behaviour of 
cosmic fluid recently reported by several authors. 

Our point of view is, in some sense, conservative since we want to investigate if and when 
such a behaviour can be achieved either using fundamental effective theories of gravity 
or by correcting standard perfect fluid cosmology as in the case of Van der Waals. 
The underlining philosophy is related to the fact that quintessence is by now become a 
sort of paradigm (like inflation!) so that every kind of scalar field (and potential) could 
be safely introduced into the game also without any fundamental physics motivation. 
Such a way of thinking is disturbing since we cannot claim to solve a new problem by 
introducing new fundamental ingredients. 

With this premise in mind, we have analyzed three schemes essentially based on "already 
known" aspects of fundamental physics. 

The first is just a slight correction: what happens if, instead of perfect fluids we adopt 

Van der Waals fluids as source in Einstein-Friedmann cosmological equations? 

This means to take into account a more realistic description of the matter content of the 

universe, capable, on the other hand, of implementing dynamics of cosmic flow where 

phase transitions are enclosed. Van der Waals fluid equation can be read as a second 

order correction (in terms of matter-energy density) of perfect fluid equation. 

It is interesting to see that for wide ranges of parameters accelerated behaviours of scale 

factor can be achieved. Besides, we get interesting matches with observational data, in 

particular the Hubble-Snela diagram and the age of the universe. 

This is not a full self-consistent model but it is a useful indication of the fact that obser- 
vations can be explained using more realistic cosmological fluid sources. 
The curvature quintessence approach is completely different. Our starting point is the 
fact that effective gravitational lagrangians turn out to be extremely useful to study the 
quantum behaviour of gravity. Such a formalism is based on the introduction of higher- 
order corrections of curvature invariants in the Einstein-Hilbert gravitational action in 
order to get renormalization, at least, at a finite number of loops. 

This scheme has given rise to several interesting inflationary models (e.g. Starobinsky 
one), in which shortcomings of standard cosmological model were cured in the framework 
of such geometrical corrections. We asked for a similar approach at today very low en- 
ergies. In other words, if remnants of primordial quantum gravity effects have survived 
until now, they could be responsible of the observed unclustered dark energy and explain 
the accelerated behaviour. We worked out some fourth-order gravity cosmological models 
and found that also in this case quintessential issues can be matched. The comparison 
with observation give suitable values of H and of the age of universe. 
The last scheme investigated is a cosmological model with torsion. Introducing spin into 
dynamic as further source of gravitational field equation is the main goal of theories with 
torsion. We wonder whether some kind of torsion field, among the several ones existing, 
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can work in order to give rise to accelerated behaviours of cosmic fluid. Our results 
confirms such hypothesis and match extremely well with observational data as discussed 
above. 

In conclusion we have proposed alternative but physically founded approaches to quintessence 
in order to see if such a paradigm can be recovered in the framework of fundamental 
physics. 

Obviously the presented models have to be improved in order to match all the quintessen- 
tial issues as the coincidence problem, the low value of (dark energy) cosmological con- 
stant and the compatibility with large scale structure. 
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